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Abstract 

Properties of the self-adjusted Monte Carlo algorithm applied to 2d Ising ferromag- 
net are studied numerically. The endogenous feedback form expressed in terms of 
the instant running averages is suggested in order to generate a biased random walk 
of the temperature that converges to criticality without an external tuning. The 
robustness of a stationary regime with respect to partial accessibility of the infor- 
mation is demonstrated. Several statistical and scaling aspects have been identified 
which allow to establish an alternative spin lattice model of the financial market. It 
turns out that our model alike model suggested by S. Bornholdt, Int. J. Mod. Phys. C 
12 (2001) 667, may be described by Levy-type stationary distribution of feedback 
variations with unique exponent a\ ~ 3.3. However, the differences reflected by 
Hurst exponents suggest that resemblances between the studied models seem to be 
nontrivial. 
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1 Introduction 



During the past decades, the financial markets around the world have become 
more and more interconnected. The financial globalization has changed the or- 
ganization structure of the stock markets by creating new risks and challenges 
for market participants. Dramatic expansion of the cross-border financial flows 
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as well as domestic flows within the countries due to very rapid increase of 
telecommunication and computer-based products have emerged. As a conse- 
quence, the world financial markets are increasingly efficient today than ever 
before. 

The integration and globalization process opens discussion about the activity 
of the central economic institutions that focus attention to financial mar- 
ket monitoring. The subsequent regulatory legislative yields asymptotically 
to dynamic balance between coalitions based on agreements of the cooperat- 
ing companies and on the other hand competitors that enhance a diversity 
of the price trends. The feedback signals of different reliability are monitored 
and daily evaluated by central economic institutions. Consequently, this af- 
fects correlated decisions of the market participants. An expected effect of the 
global self-regulation is a portfolio efficiency. 

Recently, there have been many attempts to gain a control over the dynamics 
of complex systems. The approaches based on the principles of feedback and 
adaptivity emerges as a perspective branch (Jll). The importance of the positive 
feedback mechanism in the economic context has been studied in Ref. §2). An 
approach based on the feedback control of spatially extended optical system 
has been employed in the Ref. (0). The feedback mechanism is also exploited 
in Monte Carlo (MC) algorithm introduced by us recently in Ref. (J4J). The 
algorithm is formulated as a random walk of temperature variable biased by 
feedback that mixes current and past stochastic signals, originating at the 
platform of extended statistical system. The time integrated signals are used 
to build up the actions shifting an instant temperature towards the critical 
temperature value. Due to limited memory depth, the convergence process 
is necessarily accompanied by uncertainty. Since the feedback treats a partial 
information only, many parallels with risky decisions within the business world 
can be found. Therefore, the MC feedback model may be used for simulation 
of decision-making in the stochastic environment. 

The plan of the paper is as follows. In the next section the self-adjusted MC 
algorithm is reintroduced briefly. Its relationship to non-equilibrium spin up- 
dates is discussed in Sec. 3. The Sec. 4 deals with the statistics of spin clusters 
near the critical conditions and the self-adjusted algorithm is related to cellular 
automaton models of self-organized criticality (SOC). In Sec. 5, the applica- 
tions of mentioned MC dynamics to models of financial markets is presented. 



2 The self-adjusted random walk near criticality 

Before proceeding to the introduction of financial model, we recall the salient 
points of the self-adjusted MC model (0) that is based on the random walk of 
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the temperature variable T t defined by the one-step recurrence 



r t+1 = T t + £ f Asign(F t ), (1) 

where A parameterizes a maximum step length \T t +\ —T t \, randomized by £t 
that is uniformly distributed within the (—1, 1) range. The subscript t accounts 
for MC steps per N degrees of freedom. The feedback regulatory response 

_ (E*) t - 3(E) t (E 2 ) t + 2(E)* (E 2 ) t - (E) 2 
1 T t A N TfN ' 1 ' 

mixes an influence of bias and randomness involved in the energy series {E t } 
that represents outcome of some extended statistical system. Eq.(2) approxi- 
mates a temperature derivative of the specific heat ^ [( (E 2 ) eq — (E)l q ) / (T 2 N)] 
at the canonical equilibrium (eq) ensemble defined by the constant tempera- 
ture T ~ Tt (H). In Eq.(2) the canonical averages are approximated by instan- 
taneous running averages of the energy cumulants (E p ) t , p = 1,2 calculated 
with the assistance of linear filter (Jy) as it follows 

(E*>)t = (l-v)(E p )t-i + vE! ■ (3) 



Here the plasticity parameter rj e (0, 1) reweights an energy samples E t ac- 
cording to delay. For rj — > 0, running averages vary slowly and long-memory 
processes dominates, whereas in the limit 77 — > 1 — the memorizing becomes 
faster but less efficient. 

The feedback is applied to maintain the critical regime of the 2d Ising spin 
micromodel of ferromagnet. This is defined for % = 1,2, ... ,N = L 2 sites of 
L x L square lattice characterized by the exchange energy E = —J2 nn S%Sj, 
Si G {— |, |}, where exchange coupling constant rescales an instant tempera- 
ture T t . The nearest neighbor (nn) interaction is considered. This interaction 
picture is supplemented by the periodic boundary conditions in all directions. 
The standard Metropolis single-spin-flip algorithm (H; 0) where two states are 
linked by the flip acceptance probability p acc = min{l, exp (—8Ei/T t ) } de- 
pends on the energy change SEi at ith site. The technical issue of the proposed 
implementation is that T t is hold constant during N = L 2 spin updates (i.e. 1 
MC step per sample). The example which illustrates T t dynamics is shown in 
Fig.l. As follows from Eq.(3), a filtered information serves to provide T t+1 near 
criticality. In compliance with a previous facts, the feedback role in dynam- 
ics can be summarized as follows: a) it receives series of energy {E t } values 
from extended statistical system (or only sufficient part of the series is served 
for this purpose); b) it treats the series via the linear filtering that generates 
running averages {E p ) t ; c) it transmits the regulatory signals sign(i^) via the 
temperature given by Eq.(l). 



3 



E-h 0.6 - 




150000 300000 

t 



450000 



Fig. 1. Temporal dependence of self-adjusted MC model for transient T t dynamics 
obtained for L = 90 and parameters r\ = 0.01, A = 10 -5 . Two paths [starting from 
Tt=o = 0.4, and 0.8] tend to the unique attractor T C (L) (asymptotic regime), as the 
variable t increases. 

After the transient time the thermal noise overcomes a regular bias involved 
in stochastic F t and system with feedback maintain a stationary regime. Due 
to fluctuations of T t a non-equilibrium distribution is generated. Therefore, 
the expectation mean values, worth our interest are given by 



T p ) 



m - n 



E T l> 



p=l,2. 



(4) 



-»■<> 



where n\ 3> uq. If the number of excluded steps hq is larger than the transient 
time, the average (T) estimates critical point of the extended system. Thus, 
the observation time is defined by r b s = n\ — n . In our simulations we used 



n 



io 5 , m ~ io 7 . 



3 The finite-size effects of (T) 



There are many examples of the simple non-equilibrium spin-flip stochastic 
models (0; 0) giving rise to the various dynamical phases. To simulate nearly 
equilibrium systems, the specific non-equilibrium approaches have been sug- 
gested. For instance, the Ref. (fiof ) exploits the controlled growth of percolating 
spin cluster. 

The postulation of artificial dynamics clearly places novel time scales into 
original equilibrium problem. In our case the combination of two parameters 
7] and A affect the distance between a stationary non-equilibrium statistics 
and canonical equilibrium. The basic properties of stationary non-equilibrium 
regime can be summarized as follows: 
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1. The stationary regime with a ferro-paramagnetic order and positive (T) 
is stabilized by a sufficiently slow dynamics limited to A < A tT (jj), where 
A tr (ry) is a threshold value. The irreversible drop to negative (T) (antifer- 
romagnetic order) occurs at A = A tr . 

Similar transition can be found for the coupled map lattice where anti- 
ferromagnetic order is generated spontaneously by ferro magn etic couplings 
accompanied by sufficiently fast synchronous dynamics ((ill). For L = 10, 
T) = 10" 3 we estimated A tr ~ 2.2 x 1CT 3 . 

2. The canonical equilibrium limit (T) — ► T C (L), where T C (L) is the pseu- 
docritical temperature of specific heat, is reached for the sufficiently slow 
dynamics 1 ^> A > 0, when the feedback memory is deep (1 ^> rj) and 
parameters admit the scale separation ?; > A. In addition, the stationary 
averages have to be accumulated for r obs 3> I/77. 

3. For stationary regime (see point 1) the nonzero variance ((T— {T)) 2 ) gives 
rise to the broad energy distribution in comparison to Boltzmann one. 

The question of practical relevance is whether true T c can be attained by non- 
equilibrium dynamics. The result of simulation realized for rj = 10~ 5 , A = 10~ 9 
is depicted in Fig. 2. The standard finite-size dependence of T C (L) ~ (T) = T c + 
b/L (0) (where b denotes the thermal coefficient) has been obtained from the 
fit. Without additional tuning the simulation provides estimate T r ~ 0.5679 
of the equilibrium exact value T c ex ~ [21n(l + v^)] -1 - 0.56729 {2) . 
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Fig. 2. The non-equilibrium dynamics confirms 1/L finite-size dependence of 
pseudo-critical temperature. True critical value T c is estimated from the fit 
(T) =T C + b/L for L = 10, 20, . . . , 90. 

As one can see in Fig. 3a follows that substantial finite-size corrections occur 
for a short-time memory. In that case, the temperature mean value can be 
interpolated by 



(T) = T c + - - Arf , A < T) < Vc (L) 



(5) 
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Fig. 3. a) The 77-dependence of (T) calculated for L = 5, 6, 7, 8, 9, ... , 70, 80, 90 
lattices. The mean temperature decreases with 77 for 77 < rj c (L). The upper bound 
is chosen to guarantee that (T) calculated for different lattices merge to the unique 
rf asymptotics. b) The robustness (or failure tolerance) of averages with respect to 
feedback that accumulates partial statistics for the rectangular segment [<f>L] x L, 
0<({><1 of LxL lattice. 

where cut-off r) c (L) is roughly interpolatated by r) c (L) = rj^ — 0.343/L with 
7/00 = 0.468. From the fit carried out within the range 77 e (0, ^) we have 
determined the exponent q ~ 1.5. Assuming that the lowest mean value (T) ~ 
^min ~ T 2 (L), where 77^ is achieved, for L — > 00 we may write A = (T c — 
^min) I via- The presence of the temperature lower bound is associated with the 
condition that peak of specific heat is formed. This implies the requirement of 
a residual acceptance of single-spin-flips at very low temperatures. The bound 
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Fig. 4. a) The pdf's of spin cluster for L = 10, . . . , 100. The inset shows a lattice 
size dependence of the peak position s max (L) oc L D . The simulation is carried out 
for parameters r\ = 10~ 2 ,A = 10~ 5 . b) A log-log plot of spin cluster size pdf 
considered for different 77. The exponent r obtained from the fit P(s, L) oc s~ T 
within the range s < 30. The differences in r G (1.47,2.21) are connected with 
changes (T) (see Fig. 3). 

can be explained in the frame of two- level specific heat model denoted as C 2 . 
The model takes into account only single branch of ground state —L 2 /2 (since 
the ergodicity is broken, the ground state is counted once) and L 2 -fold single- 
spin excitation of energy 2 — L 2 /2. The condition for extreme dC^/dT^^ — 
yields equation (T 2 — 1) exp(2/T 2 ) + L 2 (l + T 2 ) = 0. The last relation predicts 
a monotonous decrease of T 2 with L and solution can be simply approximated 
by T 2 ~ f° r infinite L asymptotics. To estimate T 2 we have carried out 
simulation for L = 90 and i] e (0.45, 0.55) that gives (T) from (0.25, 0.3) 
range. The mean value (T) is approximated by T 2 = 0.212 ~ = 0.222. 

The assumption about the energy gap corresponding to two lowest levels is 
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consistent with the presence of the crossing over memory depth l/rj c ~ 2 (two 
MC steps). To maintain a stationary regime the substantial memory reduction 
for r] > 7] c have to be compensated by the enlarged (T) value. We observed 
that for 7] ~ 7] c the ergodic time r e between the plus and minus ma gne tization 
stages diverges as r e oc L z , where z ~ 2 is the dynamical exponent (|l3[). If r b s 
is smaller than r e , the broken symmetry can be detected. Clearly, for rj > i] c 
the "jaggy" (T)(rj) dependence indicates the occupancy of only few lowest 
energy levels. 



4 The relationship to SOC and percolation models 



The concept of SOC proposed by Bak (14) represents the unifying theoretical 
framework of systems that drive themselves spontaneously to critical regime. 
The sand pile, forest-fire (fl5h and game-of-life (flit ) are well-known examples 
of cellular automaton SOC models. Their stationary regime can be character- 
ized by invariant probability distribution functions (pdf 's) of spatial or tem- 
poral measures of the dissipative events called avalanches. The corresponding 
pdf's include a power-law intervals with the specific non-equilibrium critical 
exponents. 

According to Ref. (Jl7f ). the stationary SOC regime implies the operation of 
inherent feedback mechanism. On the contrary, our model has explicitely de- 
fined global feedback coupled to the MC dynamics. Such coupling is com- 
pletely different from the standard models of SOC, but the arguments given 
in this section suggest that additional coupling anchors statistical properties 
indistinguishable from SOC models. 

In the next, we focus our attention to the problem of partial information 
regulation to answer the question, how important can be a concrete feedback 
form. We consider spatially truncated form of feedback defined by energy 
cumulants calculated for the rectangular [4>L] x L lattice segment, where G 
(0, 1) and [x] denotes an integer value of x. The results are depicted in Fig. 3b 
for L = 60. The simulation has confirmed that stationary regime is robust 
with respect to the partial information regulation represented by rj and 1 — 
terms. As one can be see in Fig. 3, both sufficiently small rj, (rj < rj c (L)), and 
1 — have similar influence on diminishing the stationary (T) value. 



As pointed out in Ref. (18), the pdf's of spin clusters at a second-order criti- 
cal point resembles pdf's of the spatial extent of avalanches. Admittedly, the 
correspondence between spins and clusters is not unique and thus numerous 
mapping schemes could be examined. We have utilized the well-known map- 
ping to bond percolation model (0) via Wolf's algorithm (Q to identify 
cluster without overturning. To interpret the cluster statistics, we suggest the 
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universal scaling in the form 



P(s,L) = L~ 




) 



(6) 



where s is the number of spins belonging to the same cluster. In Fig. 4 are 
depicted spin cluster pdf 's. As one can see, the distribution functions for clus- 
ters s <C s max (L) clearly obey the scaling given by Eq.6, where s max (L) is 
proportional to the size of a spanning cluster. It is associated with the peak 
-P(s m ax,£) in P(s,L) distribution. Using the scaling relation s max (L) oc L D , 
the fractal dimension D ~ 1.98 has been determined from the fit, see inset 
in Fig. 4a. Assuming an asymptotic form g(x) oc x~ T in the limit x — > we 
obtained P(s,L) oc L~ f5+TD s~ T that implies independence on L if (3 = t D. 
Surprisingly, the dependencies in Fig. 4b indicate that small-scale power-law 
scenario survives for full 77 G (0, 1) range. It doesn't matter how distant an 
assembly is from the thermal equilibrium except a spanning cluster that grows 
with increasing 77. 

Since the static bond percolation model pays no attention to temporal varia- 
tions, the alternative maps should be examined, for instance those taking into 
account a short-time memory of cluster evolution (21). The question whether 
differences in a spin - cluster map definition might affect a non-equilibrium 
exponent (3, introduced in Eq.6, is left for future research. 



5 The relationship to financial market models 

5.1 Statistical properties relevant for financial market modeling 

Performed simulations for pdf of F t defined as Pf{F) = ($F,F t ) confirm that 
flat tails exhibit within the full range of rj, see Fig. 5. Therefore, it is worth to 
check functional stability of aggregated regulatory responses 



The results of analysis depicted in Fig. 5b has validated the expected scaling 



for n < 10 with estimated exponent a ~ 1.3, where /(•) is some universal 
function. Since F t is bounded, the crossing-over from Levy to Gaussian scaling 



(7) 



t'=t 
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a) 




Fig. 5. a) The pdf of F t shows a high leptocurticity that survive at different choices 
of 77. b) The scaling property given by Eq.(8) and crossing over from the Levy to 
Gaussian pdf for n = 5, 10, 20. The fluctuations of r n> t, n < 20 merge to the unique 
scaling function for a ~ 1.3. 

has been verified, see Fig. 5b. More profound connection to the price dynamics 
is discussed in the following section. 

The dynamics of self- adjusted system can be analyzed by means of autocor- 
relation functions of the time lag t: 



r (p) _ ( F t' F f'+t ) ~ Cg) ( F ht) „ _ 1 9 (Q\ 

" <W>-«*?» 2 ' P " ' ( j 

The MC model yields nearly exponential decrease of C^, see Fig. 6. The power- 
law dependence of C^ 2 ^ which becomes more pronounced for larger L is iden- 
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Fig. 6. a) Semilogarithmic logarithmic of autocorrelation function Cj- at given 
interval demonstrate the standard exponential dependence and the long-time 
anti-correlation zone is also formed, b) Logarithmic plot of the autocorrelation func- 

(2) 

tion C t exhibits power-law regime described by the slope (— q ) = —1.5. The auto- 
correlation functions indicate the typical time scale of the order 0(l/n 2 ) = O(100). 

tilled. Assuming the form ~ t~ q we have determined q' ~ 1.5. 



5.2 The financial market model formulation 



On the basis of previous facts we have proposed model based on the two postu- 
lates: (i) under the conditions of new economy the world market entities (par- 
ticipants, companies) are forced to evolve towards the longer auto-correlation 
times. It means that hardly predicable - volatile price statistics gets stuck in 
the regime which imply highest degree of uncertainty. It might be classified 
as an edge of the chaos. Such interpretation is fundamentally equivalent with 
the dynamics of the continuous phase transitions and dynamics of informa- 
tion processing (jiih; (ii) the centralized institutions are established to satisfy 
the global information and legislative needs of market entities looking for the 
arbitrage opportunities. The action of participants leads to enhancement of 
market effectiveness in the sense of the postulate (i). 
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Our simulation have shown that presented spin model that involves feedback 
has much in common with Ising spin market models of the price dynamics (j24h . 
The feedback introduced in Bornholdt's model (J25h. which was later elucidated 
by Kaizoji, Bornholdt and Fujiwara (BKF), is linked to the local in time and 
global in lattice space magnetization m t = (2/L 2 ) J2i=i Si,t- Both models are 
inspired by the paradigm of minority game (j26h . Speaking in game theory 
term, the feedback promotes the continuing tournaments between competing 
ordered (ferro) and disordered (para) magnetic phases. To maintain a balance 
between ordered ferromagnetic and disordered - paramagnetic spin phase, the 
feedback mechanism including a global cooling or heating is applied. The pre- 
dicted temperature is simultaneously transmitted to all lattice sites. Clearly, 
the effect of such step has probabilistic nature, because a freedom exists be- 
tween the market entities which can uphold or abolish older or establish new 
coalitions, respectively. Subsequently, any decision is reflected by correlations 
between entities. 



Two principal situations may be observed from the view point of the proposed 
model. At the subcritical temperatures the clusters grow, which means that 
correlations among market entities become increase. A subcritical temperature 
can increase a mean cluster size until a spanning cluster come into existence. 
The ferromagnetic domain state can be understood as a market situation 
where competition is distorted due to monopoly, duopoly or trust formation 
effects. To prevent from highly correlated stocks, the feedback triggers the 
antitrust operations (|27l). The portfolio diversification is represented by the 
fragmentation of the spin clusters. 



5.3 The comparison with other spin models of financial markets 



The spin model presented in previous sections displays several interesting sta- 
tistical properties that are particularly relevant for stochastic financial market 
dynamics. One of our challenges is to examine cross-links between our feed- 
back model and BKF model. In BKF model the local spin behavior is given 
by the competition between the local and non-local interactions involved in 
local field h t) i = 2 (J2nn $j,t — ftSi,tl m *l)> where k is proportionality factor of 
feedback instant magnetization. The probability of ±| lattice spin states is 
given by 1/[1 + exp(=p2/?/ij)], where (3 is the inverse constant temperature. The 
difference m t +i — m t is interpreted as a logarithm of the price return. Within 
the BKF model the instant price is proportional to exp(Am t ) term, where A 
is related to the ratio of the interacting and fundamental traders. 

In further we assume the equivalence of exp(Am t ) and exp(Ai^) terms con- 
sidered for different models. The interpretation of exp(Ai^) as a price term 
is in accord with requirement that price increases when the competition is 
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bull (bear) market durations 



Fig. 7. The comparison of BKF model (solid lines) and model introduced in this 
paper for the parameters listed in the table 1. a) the rescaling factor £ is used to 
attain a full data agreement except the anomaly for small lags — Ft- b) the 
oscillatory deviations from the power-law has been observed. Despite of this, the 
rough consistency of cumulative volatility of returns can be attained by rescaling. 
c) the pdf 's of bull (bear) regime durations. 

weakened. The crucial items emphasizing the common aspects of models are 
listed in table 1 and illustrated in Figs. 7 and 8. 

We see that both models yield Levy distributions with the nearly common 
index a\ ~ 3.3 related to cumulative distributions with power law slope 
—2.3 = — a% + 1 that is consistent with the empirical data (29). The bowl- 
shaped anomaly at small \F t+ i — F t \ is probably due to sign (Ft) sharpness. The 



13 




Fig. 8. Comparison for BKF and our model, as in Fig. 7, for autocorrelations, a) dif- 
ference between correlations and anticorrelations as a key distinction, b) volatility 
autocorrelations with the rescaling t' = Wt to emphasize the similarity of models. 
However, contrary to the expectations in Ref. (|2a). our fits do not confirmed hy- 
pothesis about the power-law volatility autocorrelations in the spin market models, 
c) principal difference between diffusive (the curve with the slope 0.5 corresponds 
to BKF) and super-diffusive (slope 0.7) regime. 
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Table 1 

The comparison of the statistical characteristics of models of financial market for 
L = 90. The simulation results obtained for parameters n = 0.01 and A = 10 -5 . 



The BKF parameters k = 20, (3 = 2 are chosen to attain the intermittent regime. 
The positive Ft (also mt > 0) is the signature of so called bull market regime, while 
Ft < (mt < 0) is interpreted as a bear market regime. 
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volatility of log. returns 
pdf of volat. log. returns 
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diffusion of F t , m t 


\m t+ i - m t \ 
a\ = 3.3 
non-universal 
diffusion St 5 


\Ft+i — F t \ 
a.\ = 3.3 
non-universal 
super-diffusive 5t ' 7 


Fig. 7b) 
Fig.8b) 
Fig.8c) 


conditional probabilities 

7T (bull bull) 7T (bull bear) 
7T (bear) bull) 7T (bear) bear) 


0.978 0.022 
0.022 0.978 


0.974 0.026 
0.026 0.974 





remarkable feature of this comparison is a similarity of the exponents obtained 
for pdf's of the bull (bear) market durations. In addition, the similarity of tem- 
poral aspects can be seen in coincidence of stationary conditional probabilities 
of bull-bear switching listed in the Table 1. We see that F t feedback model 
acquires the exponential form typical for the positively correlated quantities, 
whereas the magnetization variations predicted by BKF are anticorrelated. 
Clearly, this difference can be explained providing energy cumulants temporal 
filtering. The simulation of BKF model reveals some gap between Hurst expo- 
nents f)30l ). Contrary to the standard diffusive St 0,5 behaviour which is typical 
for (\m t+ st — f^tl), the superdiffusive short-time regime {\F t+ st — F t \) ~ St ' 
has been identified for F t feedback model. 

5.4 The application to inter-stock correlations 

In the subsection certain diverse application of our model is presented. As 
one sees, the MC model may be used to understand the basic features of the 
mechanism of regulation of bilateral stock cross-correlations. The role of cen- 



15 



tral institutions, modeled via the feedback, is to perform large-scale eligible 
decisions that could enhance portfolio efficiency. These contribute to F t com- 
prising a summary of the stock prices. We assume that different stocks are 
distinguished by lattice position and sign of dominant trend at given time. 
And for now the correspondence between spins and stocks follows a simple 
majority rule: S{ = \ [Si = — |) when the sell (buy) market orders prevail dur- 
ing elementary time lag (1 MC step/site), respectively. The cross-correlations 
among directly regulated - four "nearest" stocks are related to the mean ex- 
change energy term (SiSj eim (i)) optimized by F t to be extremely susceptible 
to regulation. The long-range cross-correlations stem from indirect regulation 
of the portfolio composition near to criticality. 

The presence of correlations (and anticorrelations) between pairs of stocks play 
a key role in the theory of selecting of the most efficient portfolio composition 
of the financial goods (|22l). For discussed spin model the properties of the pairs 
of stocks can be characterized by the cross-correlation coefficient 

dj,t = — E SijStf , i, j E {1, 2, . . . , L 2 } (10) 

ranging from —1 to 1. Here i,j denote the pair of randomly chosen sites. In 
the simulations the temporal averages have been calculated for n av = 10 3 steps 
and accumulated in pdf's depicted in Fig.9. We see that the pair correlations 
are dominantly positive except a negative tail that agrees qualitatively with 
the empirical studies (J22f). For r\ larger than ~ 0.1 the negative tail vanishes 
and pdf of Cj J)t is shifted towards 1. Clearly, such "coherent" regime is typical 
for inefficient portfolio structure. 



6 Conclusion 

In conclusion, we have presented application of adaptivity concept to critical 
phenomena. Despite of the principal difficulties, insufficient understanding of 
the systematic errors which arise proceeding from non-equilibrium to equi- 
librium, we expect that work might inspire numerous applications to other 
inverse statistical problems where constrains are written in terms of running 
averages. 

We have demonstrated that MC model can be reinterpreted in macro-economic 
terms. The simulations indicated spontaneous emergence of the power-law 
distribution for feedback regulatory signals which is universal feature of the 
stationary non-equilibrium regime. Surprisingly, we found that F t+ i — F t statis- 

(2) 1 

tics contradicts to our previous expectations based on the finding C t ~ t~ q , 
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Fig. 9. Cross-correlations distribution obtained for different r\ [see Eq.(10)]. The 
anticorrelations appear for rj = 0.01. In that case the distribution can be approx- 
imated by double-Gaussian P{c) = J2k=o 1 w k ex P [ ~~ ( c ~ a k )V (20^) ] including 
the weights wq = 1.432, w\ = 1.951 and dispersions ctq = 0.146, o\ = 0.096 and 
a = 0.337, ai = 0.493. 

which does not imply power-law fall-off of the autocorrelations of \F t+ \ — F t \. 
The influence of the memory depth related to parameter rj has also been in- 
vestigated extensively. The critical point related to the most efficient portfolio 
structure, is an abstract of ideal situation when rj — > + , which is unattainable 
in a real finite-memory conditions. Otherwise, we have found that the portfo- 
lio efficiency is weakened via short-time memory for rj > 0.5. The comparative 
study of two spin models have confirmed that presented MC dynamics based 
on the specific heat feedback offers an alternative description of financial mar- 
kets. Despite of outstanding similarities, the comparison of exponents indicate 
non-trivial relationship between the universality classes of studied models. 
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